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Abstract. - We propose a model with a quantized degree of freedom to study the heat 
transport in quasi-one dimensional system. Our simulations reveal three distinct temperature 
regimes. In particular, the intermediate regime is characterized by heat conductivity with 
a temperature exponent 7 much greater than 1/2 that was generally found in systems with 
point-like particles. A dynamical investigation indicates the occurrence of non-equipartition 
behavior in this regime. Moreover, the corresponding Poincare section also shows remarkably 
characteristic patterns, completely different from the cases of point-like particles. 



The problem of heat conduction in low-dimensional systems has been an interesting sub- 
ject for a long time due to the potential applications of nanostructured materials and the 
importance in understanding the physics of low-dimensional systems (see review [1] and the 
references therein). A long standing issue is the exploration of the necessary condition for 
a dynamical system to obey the Fourier law [2-6]. For this purpose, a class of "Billiard 
gas channel" has been proposed recently [7-10] to investigate the connection between macro- 
scopic heat transport and microscopic dynamics. In these models, usually a series of scatterers 
are placed periodically along two parallel walls, and point-like particles without interaction 
frequently undergo specular reflections on the boundary formed by the scatterers and the 
walls. By tracing the motion of the particles, the energy transport in a periodic structure 
is investigated. Despite intensive simulations in recent years, the sufficient condition for a 
dynamic system ensuring the Fourier law is still not clear. For example, it was ever found 
that the system with positive Lyapounov exponent [7] obeys the Fourier law, while the same 
result was obtained in some linear mixing systems with zero Lyapounov exponent [8-10]. Re- 
cently, a modified Lorentz model where particles collide with fixed freely-rotating scatterers 
was suggested to investigate the heat transport in local thermal equilibrium [11,12]. Normal 
transport is obtained in such system, indicating that the transport coefficients are finite in 
the thermodynamic limit. 
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Although these simple models provide a framework in which heat transport emerges, they 
assume point-like particles (heat carriers), and fail to provide a clear understanding of tem- 
perature effects on heat transport. In fact, along with the dramatic achievements in nanotech- 
nology, the control of heat transport for nanoscopic devices becomes increasingly important. 
Thus, it is crucial to understand the temperature dependence of heat conductivity in low 
dimensional system. Since the billiard model captures the underlying dynamics of the system, 
it will be appropriate for us to study the heat conductivity from the microscopic dynamic 
point of view. 

In this paper we propose a simple model to study the heat transport, which demonstrates 
different regimes in the heat conductivity. We consider a quasi-one-dimensional billiard gas 
channel in which the heat carriers possess an additional (intrinsic) quantized degree of free- 
dom (e.g., rotational velocity) beyond two translational ones. In general, this kind of degree 
of freedom can come from the composite nature of the carriers, or from the confinement in 
a third dimension. We explore the scattering of the particles with rough surfaces (modeled 
by fixed half disks). Here we focus on inelastic events and take no account of the energy or 
the momentum exchange between the particles and the scatterers. Namely, a collision with 
the rough surface involves the occurrence of energy transfer between the translational and the 
rotational degrees of freedom. Our simulations show that there are three regimes according 
to the temperature exponent 7 in the heat conductivity k ^ T"^. The transition from high 
temperature regime to intermediate temperature regime occurs close to a characteristic tem- 
perature Tc: below which the translational velocity distribution deviates from the Boltzmann 
distribution while energy is distributed unequally among the three degrees of freedom. In 
addition, by plotting the Poincare section we find the clear emergence of strange attractor 
accompanying the frozen of rotation degree of freedom at low temperatures. 



Fig. 1 - Schematic of the collision of a particle of radius r with a scatterer in a unit cell of the channel. 
The surface of the semicircle scatterers are plotted by dashed lines, near which the solid lines are the 
surrounding tracks of the center of particles when colliding with the scatterers. In the calculations, 
L = 3.6a, W = 1.0392a, R = 0.89a and r = O.Ola. 

In our model, the channel in x direction consists of two parallel walls of distance W (Fig- 
ure ^ shows the unit cell of the channel). A series of semicircle scatterers with radius R are 
placed periodically along two horizontal walls. Each cell has two semicircles, one on the bot- 
tom, the other on the top. The heat baths at two terminals of the channel are modelled by 
stochastic kernels of Gaussian type [9] , 
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which determines the translational velocity Vx and Vy of the particles in the channel after 
colliding with the walls between the channel and the heat baths. 
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The particle with radius r has a rotational velocity uj in addition to the conventional 
velocities Vx and Vy which describe the particle's translational movement. Since each unit 
cell contributes only one heat carrier whose radius is small compared with the scale of the 
channel, the interaction between the particles is not taken into account for simplicity. Under 



the assumption of v'^ 



, the conservations of energy and angular momentum read 
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where the primed and unprimed ones represent the quantities after and before a collision; 
the v\\ and v± refer to the components tangential and perpendicular to the scatterer surface, 
respectively. The moment of inertia around an axis through the center of a particle of mass 
m is / = amr^. The similar collision rule has been appeared in [11,12]. Note that in our 
model a ranging from to 1 measures the moment of inertia, which corresponds to the mass 
distribution of the disk particle, i.e., a = for point particles, 1/2 for uniform distributed 
mass. Eqs. Q and Q lead to the following solution: 
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In addition to the specular reflections with the parallel walls, the collision rules of semicir- 
cles are described as above equation. Both the translational and the rotational velocities are 
changed after a collision owing to the friction force against the total tangential velocity of the 
contact point at the surface of the particle (ball). 

The rotational degree of freedom naturally introduces the quantization condition luji = 
Ih, I = 0, ±1, ±2, • ■ • which is well known for angular momentum in quantum mechanics. 
Clearly, the moment of inertia determines the spacing of the quantized angular velocity Aw = 
h/I = h/amr^, and the energy spacing h^/amr"^ between the ground state and the first 
excited state defines a characteristic temperature T^, = / amr^kB- Actually, the rotational 
part can be regarded as an intrinsic degree of freedom, then our model itself manifests internal 
excitations. If the energy for such excitation is much less than ksTc, the internal degree of 
freedom is frozen. In our simulations, the particle's ui takes an allowed discrete value near the 
w' given by eq. . Under this assumption, the conservation of angular momentum is broken 
down, for the tangential component of velocity v\\ is the same as , while the perpendicular 
component v~i is determined by the energy conservation. In the simulation process, we set the 
parameter a, the particle mass m, Boltzmann constant kB and Planck's constant fi as unity 
for convenience. 

As the quantized degree of freedom is assumed to play the role in the channel, the rotational 
velocity then does not change when a particle in the channel is rebounded from the walls of 
two heat baths. The heat flux is calculated by the exchange of energy carried by the particles 
via the collisions with a heat bath within total time t m spent for M such collisions [8] , 
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where {AE)j = {Ein — Eout)j is the energy exchange at the jih collision with a heat bath. 
Note that above heat flux is for a single particle. In the channel with N cells, the heat 
conductivity is found to be k = JN'^L/AT. 
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Fig. 2 - (a) Temperature dependence of heat conductivity, (b) The ratio of average rotational 
energy over average translational energy at different temperatures. The calculations are made for 
a = 0.1,0.2, and 0.4. 

Figure |2Ia) sho-ws the heat conductivity-temperature curves for a channel with = 32 
cells. Here = 5000/a (in units of /mkBO?). Because of non-equipartition of energy in the 
system, "which "we "will discuss next, one can not define a local temperature from microscopic 
point of vie"w. For simplicity, "we assume an average value of temperature T = (Tl + Tb)/2 
under small temperature gradient, -where Tl and Tr are the temperatures of t"wo heat baths 
respectively. The temperature difference AT/T is fixed at 0.2. To obtain the heat conductivity 
correctly, "we investigated heat flo"w J in Eq. jSJ as a function of temperature gradient in a "wide 
range of 5 orders of magnitude, as "was done in ref. [13]. We have verified the linear relation of 
J ~ VT for temperatures and gradient explored in fig. El One can clearly see that the curves 
have larger slope at temperatures near and less than Tc. Our calculations exhibit that the 
maximum value of the slope reaches 1.422 for the case of a = 0.4. 

To illustrate this characteristic, "we investigate the average translational energy {Etji and 
the rotational energy {Er)i in a unit cell 

= - = t^r. (6) 

"where Et = ^mv^ and Er — i represents the ith cell and tij is the time spent "within 

the cell for the jth visit. In fig. IHb) "we plot the ratio of {Et)/{Er), "where (Et) and (Er) 
are the translational and the rotational energies averaging over N cells, respectively. Clearly, 
there are three regimes "with different characteristics. 

1. Lo"w temperature limit: Our simulations exhibit k ~ T'^ "with 7 1/2 for various values 
of a at lo"w temperatures (sho"wn in fig. I2a)). This coincides "with the result in previous 
literature [8]. We also calculate the temperature exponent 7 for the model of [7] and [14], 
respectively, obtaining the result 7 ~ 1/2 again. In the follo"wing we make an estimation of 
the heat conductivity k. In fact, at lo"w temperatures T ^ Tc, the heat carriers in our model 
can be regarded as point-like particles because their thermal energy ksT is much smaller 
than the first excitation energy. In this case, the rotational motions are frozen. The heat 
conductivity is estimated as k ^ X{\v\) "where A is the mean free path and (\v\) describes the 
average absolute velocity [15]. For ideal gas it follo"ws that {\v\) ^ yJ^k^Tpm. Consequently, 
"we have k ~ T^l"^ when A is treated as a constant. 
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2. Intermediate temperature regime: The abrupt increase for the curves in fig.|2Ia) in the 
intermediate temperature regime is related to the presence of quantized rotational degree of 
freedom. It seems that the less moment of inertia the particle possesses, the further the abrupt 
increasing regime is away from its characteristic temperature. As is shown in fig.EIb), the ratio 
(Et) / {Er) rises rapidly from zero with increasing temperature T till Tc, then asymptotically 
to 1/2. This indicates that the total energy of the system is not equally partitioned among 
the three degrees of freedom except in high temperature limit. As a result, the temperature 
dependence of the average velocity {\v\) is significantly affected by the ratio of the translational 
kinetic energy to the total energy, which seems to be responsible for the abrupt increase of 
heat conductivity since k is proportional to {\v\). The non-equipartition of energy, which 
has been studied in granular mixtures [16-18], has now been observed in our model due to 
non-conservation of angular momentum. 

3. High temperature limit: In high temperature limit, the saturation of the value of 
{Et) / (Er) is observed. In this case, the quantum effect are almost overwhelmed by the thermal 
fluctuations and thus the macroscopic system is in thermodynamic equilibrium. The curves in 
fig.[5Ja) show 7 goes back to 1/2 again, however, the underlying mechanism is quite different. 




Fig. 3 - Probability distribution at temperatures T/Tc = 0.8 (o), 2 (•), 4 (o), and 10 (★) for a = 0.4 
and N = S. (a) The translational velocity distribution. Solid lines are the Boltzmann distributions 
at the corresponding temperatures, (b) The rotational velocity distribution. The T* is obtained by 
the best fits (solid lines) of Gaussian distribution for numerical data. 

In the non-equilibrium regime, the traditional definition of temperature can not well de- 
scribe the features of the system. To understand the characteristics of these regimes further, 
we measure the translational velocity distribution P{v) and the rotational velocity distribu- 
tion P{u}) of the particles. The simulation is performed in a equilibrium system of = Tr, 
where 10^ particles with initial Gaussian distribution of velocity at a certain temperature T 
are put into the channel to make collisions with the semicircle rough surfaces, while under- 
going specular reflections with the straight walls. After a long relaxation time, the system 
reaches a steady state. As shown in fig. ^a), solid lines are the Boltzmann distributions in 
equilibrium at the corresponding temperatures. Better agreements appear when T/Tc > 10, 
which indicates that the particles can be considered as the classical ideal gas in equilibrium 
in the high temperature limit despite the existence of quantized rotational velocity. In com- 
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parison to the velocity distribution at relatively low temperatures, the deviations from the 
Boltzmann distribution are evident, leading to a larger average velocity. Figure |2Ib) shows 
the rotational velocity distributions P(Jjj), which are typical Gaussian distributions. T* is the 
best fits of the curves. As expected, it is lower than the given temperature. The result of 
T* /T 1 with increasing temperature till the high temperature limit shows that the system 
is asymptotically close to equilibrium state. 
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Fig. 4 - (Color online) Poincare surface of section projected to s — v\\ plane, by following the first 
10** bounces with zero initial angular velocity. The particle of velocity v and a — 0.4 starts from the 
center of the left boundary of the cell with an incident angle of 0.9. (a) v = 100, and the allowed 
o) = 0; (b) 11 = 200, and the allowed Cj = -25000, 0, and 25000; (c) v = 1000, and more values of Cj 
are allowed. Additionally, the Poincare surfaces of section for ub — —25000, and 25000 at u = 200 
are displayed on the second row, respectively. 

To capture the physics of our model, we investigate the Poincare surface of section (SOS) 
which has been employed to characterize the microscopic dynamics for the motion of a point- 
like particle colliding with the boundary [19]. It is thus natural to investigate the dynamical 
characteristics of the system in the presence of quantized degree of freedom. Using our previous 
strategy [14], we treat one unit cell of the channel and close the two ends by straight frictionless 
walls. Particles with certain initial velocity are injected into the cell, then they collide with 
the boundary. The SOS in our model is defined as (s, v\\ , uj) where s is the distance to 
the starting site along the boundary, t>|| the tangential component of translational velocities 
and Lo the rotational velocity of the particle. By projecting the points in phase space to 
the s — ?J|| plane, we plot the Poincare surface of section (s, tJ||) where different values of 
rotational velocity uj are displayed in different colors. In our simulations, the initial incident 
angle brings about minor effects on the SOS. Figure 01 (a), (b), and (c) show three distinct 
cases of w = 100, 200, and 1000, respectively. As a comparison, we investigate the SOS for 
particles with a classical rotational degree of freedom and find it is filled with widely dispersed 
points. For the case oi v — 100 (shown in figs-^ja)), the rotational velocity is frozen, thus 
the points are confined in the plane & — 0. Comparing the pattern with that of classical 
situation, we note that many areas are eclipsed. The remained pattern implies the existence 
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of a fractal attractor [20,21]. This is completely in agreement with what occurred in Bohr 
model of atoms. As is known that the admissible orbits of an electron moving in the Coulomb 
field of nucleus occupy the whole phase space. The Bohr-Sommerfeld quantization condition 
excludes many of the classically admissible orbits, which results in the emergence of eclipse 
in phase space. When v — 200, the rotational velocity is excited and the allowed oj takes the 
values of —25000, 0, and 25000. The corresponding SOS are shown in the bottom panels of 
fig. 21 The points are dispersed in the three layers and form characteristic patterns in each 
layer. For the case oi v = 1000, more excited states of higher rotational energy are allowed, 
as shown in fig.E^c), and the points at different uj are randomly dispersed. 

In summary, we have proposed a method to introduce the quantized degree of freedom 
for the first time in the study of heat transport in the billiard channel. We have demon- 
strated three characteristic regimes for heat conductivity in the presence of quantized degree 
of freedom. Our results show that the temperature-dependent interplay between the trans- 
lational and rotational degrees of freedom plays a crucial role in heat conduction. We have 
also investigated the dynamical properties in these regimes. The non-equipartition of energy 
between the translational and rotational degrees of freedom is explicitly shown except in the 
high temperature limit, and strange attractor occurs remarkably in the low temperature limit. 
We expect that our model provides a strategy to study the effects of the quantized degree of 
freedom on heat transport. 
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